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Introduction 


Ising spin model 


Hamiltonian for a spin system 


on the lattice aZ?: o + ¥ 
ae 


i,j: labels for lattice points 
g- +1: spin variable . 


Let us consider the ferromagnetic case: Jij > 0. 
ə Spins are aligned parallel = Energy H takes lower values. 
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Ising spin model 


Hamiltonian for a spin system 


on the lattice aZ?: t + 7 + y 
= So ISh far Pal Ki ri 
ji gza 4 Yy y 


i, j: labels for lattice points 
Si = +1: spin variable . 


Let us consider the ferromagnetic case: Jij > 0. 
ə Spins are aligned parallel = Energy H takes lower values. 


What is properties of the following Gibbsian measure? 
u({Si}) = exp(-—H/T +h X S;)/Z. 
(Z = Trexp(—H/T + hX; S;): partition function) 
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Introduction 


Phase structure of the Ising model 
Free energy of the system: F(T,h) = -T ln Z = —T In Trexp(—H/T+h)-, Si). 
Rough estimate on F(T, 0): 


F(T,0)~ E-T nW. 


(W: Number of spin alignments with H({S;}) = E(:= (H))) 
ə T — œ: The second term becomes dominant, and spins are randomized. 
e T — 0: The first term becomes dominant, and spins like to be aligned. 
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Phase structure of the Ising model 
Free energy of the system: F(T,h) = -T ln Z = —T ln Trexp(—H/T+h)°, Si). 
Rough estimate on F(T, 0): 
F(T,0)~ E-TInW. 


(W: Number of spin alignments with H({S;}) = E(:= (H))) 
e T — co: The second term becomes dominant, and spins are randomized. 
ə T — 0: The first term becomes dominant, and spins like to be aligned. 


Phase structure of Ising spins: 


Te T 


For T < Te, there exists discontinuities in OF'/Oh when crossing the blue line 
(1st order PT line). 
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Introduction 


Magnetization 
ain). 
M(T,h) = a 7 (Si). 


At T < Te, M jumps as crossing h = 0 (1st order phase transition. ) 
At h = 0, naive Gibbsian measure u is not well defined 
= The system must be specified with the boundary condition at infinities: 


wT, h = 0,a) = au, (T) + (1 — a)p_(T) (0<a<1). 
with w+ (T) = WT, h> +0) 
M M 


C 


h 


F 
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2nd order phase transition 
At T = T,, there exists no discontinuities on first derivatives of F(T, h). 
What about second derivatives? = Magnetic susceptibility: 


or 


t 


As T — Te + 0, the susceptibility diverges, 


x ~ |(T — Te) /Te| 7 4 œ. 
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2nd order phase transition 
At T = T,, there exists no discontinuities on first derivatives of F(T, h). 
What about second derivatives? = Magnetic susceptibility: 


or 


As T — Te + 0, the susceptibility diverges, 
x~ |(T — T.)/T.|~7 > ov. 
Other scaling properties: 
OF 


aT? 
Mx |(T-T,)/T.|f ou 


C := -(T-T;) 


~ |(T = Te) /Te| 7% + 00, 


Scaling relation (Rushbrook identity): 
at+2b+7=2. 


(Only two of scaling exponents are independent!) 
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Scaling hypothesis 
Assume that the free energy G(T — Te, h) = F(T,h) satisfies (Widom, 1965) 
G(A*t, Ah) = Gt, h)/d. 


Taking derivatives of the both sides, 


M(t,h) = AĦIM(A t, 7A), 
x(t,h) = MTthY(AE, ATA), 
Clt,h) = d28+10(A%¢, A"R). 


Putting h = 0 and À = 1/|t|S, we get 
M(t,0) ~ jt FD /s, x(t,0) ~ |t}-@rt b/s, C(t,0) ~ |e} Gets, 
Then, a+ 26+ y = 2 because 


2s+1 
Q = 
S S 8 
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Introduction 


Block spin transformation (1) 
What justifies scaling hypothesis? = Renormalization Group 


FZ 
= 


The number of total degrees of freedom: N ++ N’. This defines the scaling factor 


Block spin transformation (Kadanoff) 


A= N'/N(< 1). 


Idea: Long range feature is difficult to be computed from microscopic theory. 


Important! 


Let's perform coarse graining, and compute correlations between averaged spins! 


Yuya Tanizaki (University of Tokyo, RIKEN) FRG & its applications May 10, 2014 @ IT Chiba 9 / 38 


Introduction 


Block spin transformation (II) 


We can generally denote the BST as 


RASP = DO TSE}, £5:})n({5:}).- 
{Si} 
T is a transition probability, which satisfies 


SE seh 


{si} 


Assumption: upr is also Gibbsian with respect to block spins s®, 
UR X exp—Hr. 


Hr(S®) = RH(S) := -ln uR. 
Free energy per unit cell: F(Hpr) = F(H)/X. 
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Introduction 


Renormalization group flow 
Operation R can be performed repeatedly on effective Hamiltonians H. 


Critical manifold 


VO ee 


Renormalized trajectory Infinitely massive 


theory 


Fixed point 


Fixed points H* = RH*=> The system shows universal /self-similar behaviors 
(Wilson). 


Important! 
RG provides a useful framework to extract and treat large-scale behaviors. | 
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Aspects of FRG 


Functional renormalization group 


FRG & its applications 


DS sles sey) General framework 


General framework of FRG 


Generating functional of connected Green functions: 


exp(WJ]) = [re exp (—S[®] + J - ©). 
infinite dimensional integration! 


Possible remedy: Construct nonperturbative relations of Green functions! 
(= Functional techniques) 


ə Dyson-Schwinger equations 
@ 2PI formalism 


ə Functional renormalization group (FRG) 


Yuya Tanizaki (University of Tokyo, RIKEN) 


FRG & its applications May 10, 2014 @ IT Chiba 13 / 38 


POSA General framework 


Wilsonian renormalization group 


Classical action/Hamiltonian: S[¢] = f dfx [ivo + mg? 4 Agdl. 
“Block spins”: ¢p = f d¢xe~!?" (x) for small momenta p. 
Field theoretic formulation of RG (Wilson, Kogut 1974): 

exp —Sa[@] := Na f II dop exp —S[¢]. 


|pI=A 


Important! 


Correlation functions (py, +++ @p,,) with low momenta |p;| < A are calculable with 
Sa instead of the microscopic action S. 
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Sonal aries 
Functional/Exact RG (1) 


Schwinger functional W; with an IR regulator Rg: 
1 
exp(W; [J )= | Does (-s Hl- 56 Reo + 5-4). 


Rx: IR regulator, which controls low-energy excitations (p? < k?). 
k-derivatives of the both sides: 


On exp(W;[J]) 


[20- 30nn (Slo) - 30o + 30) 


16 ô 
= -357 oT exp W,.[/] 


Flow equation 
16W, „ OW, 1, Wk 


rW, = Rasy — hte ayes 


2 ôJ 
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Sonal rasa 
Functional/Exact RG (II) 


The 1PI effective action Tis introduced via the Legendre trans.: 
1 
Daly] + 50° Ree = Jiel: p — Welle), 
which obeys the flow EQUATION (Wetterich 1993, Ellwanger 1994, Mortis 1994) 


OR 


=. 1 O71 Rr = 
Ont el ®] = 9ST Sra /6050 + Ry 


Properties of T: Tk > S as Rk > co, and rk >T as Rk > 0. 


Important! 


Functional implementation of “block spin transformations” keeps 
all the information of microscopic systems. 
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ORSA General framework 


Generalized flow equation of FRG 


5S;,[®]: Some function of ® with a parameter k. (IR regulator) 


k-dependent Schwinger functional 


Aspace of Noninteracting 
effective theory 
exp( (Wil ) = [De exp [— (S[® ] =r OS; [®]) + J ° ®| action 
Flow equation 


—O;,W,[J] = (On0S4[®]) 7 Interacting theory 
= exp (—W,lJ]) 3k (53k) [5/5] exp(WelJ]) >- 


C 


onsequence 
We get a (functional) differential equation instead of a (functional) integration! | 
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Aspects of FRG Optimization 


Optimization 


Choice of IR regulators ôS% is arbitrary. 


Optimization: 
Choose the “best” IR regulator, which validates 


systematic truncation of an approximation scheme. 


Optimization criterion (Litim 2000, Pawlowski 2007): 


e IR regulators ôS% make the system gapped 


by a typical energy k?/2m of the parameter k. 


ə High-energy excitations (> k?/2m) should 


decouple from the flow of FRG at the scale k. 


Aspace of Noninteracting 
effective theories k 00 
actions 


Interacting theory k > 0 


ə Choose ôS% stabilizing calculations and making it easier. 
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Aspects of FRG Examples of FRG 


Conventional approach: Wetterich equation 
At high energies, perturbation theory often works well. 
=> Original fields control physical degrees of freedom. 


IR regulator for bare propagators (~ mass term): 6S;,[®] = 16, RO’ Os. 
E à^ E^ 


EF 


— A 


Flow equation of 1PI effective action r [®] (Wetterich 1993) 


OR, 
OnRr _ 
621 ;,[®]/d@6 + Ry 


1 
OnT k [E] = 2 STr 
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Dares sinks 
FRG beyond the naive one: vertex IR regulator 


In the infrared region, collective bosonic excitations emerge quite in common. 
(e.g.) Another low-energy excitation emerges in the ®® channel 


Vertex IR regulator: 6S; = 4 gp? &, Gg) Ds. 
E} ; 


Flow equation with the vertex IR regulator (YT, PTEP2014, 023A04) 


ate] = w+ O OOA +Q + 
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Application of fermionic FRG to the BCS-BEC crossover 


Application of fermionic FRG to the BCS-BEC crossover 


FRG & its applications 


Cold atomic physics 


Ultracold fermions provides examples of strongly-correlated fermions. 
High controllability can tune effective couplings with real experiments! 


T T 
2L BEC phasea > 0) | 
T Bound state | 
E 
bt Bae) | 
T i 
5 
S 5 |/ BCS-BEC 
2 p | crossover 
uw 5 
LP L l BCS phase a < 0) | 
f No bound state 
al | ] 
, i i 
Interparticle distance F = Magnetic field [Gauss] 2000 


(Typically, T ~ 100nK, and n ~ 1011714 cm~?) 
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BCS-BEC crossover 


EFT: Two-component fermions with an attractive contact interaction. 


S= fate Fe) (2. Tae ) ER AOAO OTAG 


-0-.. 
r AS, 
ó Š, o = 
i an e 
i ] ° % 
i i © 
ka x6 2? a” 
ÅĄ& 
BCS Cooper pairs Unitary gas Dimers BEC 
< | 
0 (kras)! (Eagles 1969, Legget 1980, 


Nozieres & Schmitt-Rink 1985) 


Question 
Is it possible to treat EFT systematically to describe the BCS-BEC crossover? | 
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Application of fermionic FRG to the BCS-BEC crossover 


General strategy 


We will calculate T./ep and u/ep. 
= Critical temperature and the number density must be calculated. 


We expand the 1PI effective action in the symmetric phase: 


rey] = BFB, u) + | PG p) — Ue(p)]bp 


1 
pon © "7 J G=1(p) = Xo(p)’ 
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Application of fermionic FRG to the BCS-BEC crossover BCS side 
: 
BCS side 


Case 1 Negative scattering length (kra,)~' < —1. 
=> Fermi surface exists, and low-energy excitations are fermionic quasi-particles. 


Shanker’s RG for Fermi liquid (Shanker 1994) 


ky 
y 
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Application of fermionic FRG to the BCS-BEC crossover BCS side 


Functional implementation of Shanker’s RG 


E 4 
RG must keep low-energy fermionic excitations 
under control. 
= 58 = J, Dp RE (P) Yp with Pml 
EF < > 
k? k? s P 
RP (p) = senelo)) ( $- = eI) 6 (5 ~ Ie!) £ 
m m le 


Flow equation of the self-energy “i, and the four-point 1PI vertex r®. 


ôk +f} = a On = DE 


Ti 
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Application of fermionic FRG to the BCS-BEC crossover BCS side 


Flow of fermionic FRG: effective four-fermion interaction 


ə Particle-particle loop = RPA & BCS theory 
ə Particle-hole loop gives screening of the effective coupling at k ~ kp 


0.1 T T T T 
PP flow 
PP+PH flow 
0.08 + J 
& 0.06 F 
= 
x 
© 
ș 0.04 F 
0.02 + 
0 


4 (YT, G. Fejős, T. Hatsuda, 
km p arXiv:1310.5800) 


yp 
TBOS = ep 82-7 /2kr las] -> TBCS /2.2. (Gorkov, Melik-Barkhudarov, 1961) 
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BOSE 
Flow of fermionic FRG: self-energy 


Local approximation on self-energy: Uz (p) = ox. 


(kpa,)" 2D DA 


= 
ZS ost 
6 

0.6 b 

osh 

o ! 1 1 1 1 
0 1 2 3 4 5 8 (YT, G. Fejés, T. Hatsuda, 
kN2m arXiv:1310.5800) 


ə High energy: ox œ (effective coupling) x (number density) ~ 1/k 
@ Low energy: kak ~ 0 due to the particle-hole symmetry. 
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Bom 
Transition temperature and chemical potential in the BCS side 


(YT, G. Fejés, T. Hatsuda, arXiv:1310.5800) 


— ; j j 7 Baie 
0.45 F PP+PH a 
PP+PH+SE 


Toler 
Wee 


PP+PH+SE 
3 


-0.6 -04 02 o 2 AS 


kas) 1k ae) 


Consequence 


ə Critical temperature T../€r is significantly reduced by a factor 2.2 in 
(kras) + < —1, and the self-energy effect on it is small in this region. 


ə u(Te)/e£r is largely changed from 1 even when (kras)~' < —1. 
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Application of fermionic FRG to the BCS-BEC crossover 


BEC side 


Case 2 Positive scattering length : (kpa,)~! >> 1 
=> Low-energy excitations are one-particle excitations of composite dimers. 


atom 

Oo 
N dimer 
E E = $8 

atom 
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Application of fermionic FRG to the BCS-BEC crossover BEC side 
. 
BEC side 


Case 2 Positive scattering length : (kras)! > 1 
= Low-energy excitations are one-particle excitations of composite dimers. 


atom 
Oo 
N dimer 
© --------> 
E Es = y 
atom 


Several approaches for describing BEC of composite bosons. (Pros/Cons) 
ə Auxiliary field method 
(Easy treatment within MFA/ Fierz ambiguity in their introduction) 
ə Fermionic FRG (< We develop this method!) 
(Unbiased and unambiguous/ Nonperturbative treatment is necessary) 
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Application of fermionic FRG to the BCS-BEC crossover BEC side 


Vertex IR regulator & Flow equation 


Optimization can be satisfied with the vertex IR regulator: 


PREP) re 
s= fH [3 ee ee 
1— gR! 1— gR®(p) +3 q 455 q 43 q Ta q 


Flow equation up to fourth order (YT, PTEP2014 023A04, YT, arXiv:1402.0283): 


~< 


ôk -W = Ok = + 


Effective boson propagator in the four-point function: 


1 mas [. p? 
= (æ + z) — Re) (p) 


4m 


Tr (p) Sr 
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Begs 
Flow of fermionic FRG: self-energy 


Flow equation of the self-energy: 


OT (p+) 
x = ` 
kk (Dp) fa +12/2m + 1/2ma2 — Xz (l) 


If [Ex (p)| < 1/2mas, 


/ re (p+) 
il? + 2? /2m + 1/2ma? 
2 
J dq (8m/m?a,)np(q?/4m + =$ RO (q)) 
(27T)? ipo + a + mèas RØ (q) _ (+p)? _1ı 


2m 2ma2 


Xr(p) 


Estimate of |£; (p)|: 


DAWIDS : 5 x (V2mTa,)? x ng(k?/4m). 


2ma 


=> Our approximation is valid up to (k?/2m)/T ~ (kpas)? < 1. 
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Bega 
Critical temperature in the BEC side 


Number density: 


7 —2 
n = / ip? + p?/2m + 1/2ma2 — Xo(p) 
QmT,)3/2 [T 


Critical temperature and chemical potential: 
T./ép =0.218,  u/ep = —1/(kpa,)?. 
= Transition temperature of BEC. 


Consequence 


FRG with vertex regulator provides a nonperturbative description of many-body 
composite particles. 
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BCS-BEC crossover 
fermionic FRG for the BCS-BEC crossover 


We discuss the whole region of the BCS-BEC crossover with fermionic FRG. 
= Combine two different formalisms appropriate for BCS and BEC sides. 


Minimal set of the flow equation for Xi, and TO vt, arXiv:1402.0283) 


ôk -W = s ae 
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BGS BEE mesic! 
Flow of fermionic FRG with multiple regulators 


Flow of four-point vertex: 

Important property: fermions decouple from RG flow at the low energy region. 
ə In BCS side, fermions decouples due to Matsubara freq. (k?2/2m < TT). 
ə In BEC side, fermions decouples due to binding E. (k?/2m < 1/2ma?). 


Approximation on the flow of the four-point vertex at low energy: 


Or 


K 


Flow of self-energy: 
At a low-energy region, the above approx. gives 


On — > = + 
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Qualitative behaviors of the BCS-BEC crossover from f-FRG 
Approximations on the flow equation have physical interpretations. 


Four-point vertex: Particle-particle RPA. The Thouless criterion 
1/T® (p = 0) = 0 gives 


~=-2[" Jones ae 5(e— u) J 


2(e — u) 2e 


=> BCS gap equation at T = T,. 
Number density: n = —2 f 1/(G~+ — X). 


(T) a [D Ata, mA 
n= -2 : co) -5, | In ji+ COR 


= Pairing fluctuations are taken into account. (Nozieres, Schmitt-Rink, 1985) 


Consequence 
We established the fermionic FRG which describes the BCS-BEC crossover. | 
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Summary 


Ae 


FRG & its applications 


Summary 


Summary 


ə RG provides a useful framework to extract and treat large-scale behaviors. 


ə Functional implementation of coarse graining provides systematic treatment 
of field theories. 


e@ Fermionic FRG is a promising formalism for interacting fermions. 
=> Separation of energy scales can be realized by optimization. 
= Very flexible form for various approximation schemes. 


e@ Fermionic FRG is applied to the BCS-BEC crossover. 
=> BCS side: GMB correction + the shift of Fermi energy from u. 
= BEC side: BEC without explicit bosonic fields. 
= whole region: Crossover physics is successfully described at the 
quantitative level with a minimal setup on f-FRG. 
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